Vertex renormalization in DC conductivity of doped chiral graphene 
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The remarkable transport properties of graphene follow not only from the the Dirac-like energy 
dispersion, but also from the chiral nature of its excitations, which makes unclear the limits of 
applicability of the standard semiclassical Boltzmann approach. In this paper we provide a quantum 
derivation of the transport scattering time in graphene in the case of electron-phonon interaction. 
By using the Kubo formalism, we compute explicitly the vertex corrections to the DC conductivity 
by retaining the full chiral matrix structure of graphene. We show that at least in the regime of 
large chemical potential the Boltzmann picture is justified. This result is also robust against a small 
sublattice inequivalence, which partly spoils the role of chirality and leads to a doping dependence 
of the resistivity that can be relevant to recent transport experiments in doped graphene samples. 

PACS numbers: 72.10.Di,63.20.kd,81.05.Uw 
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I. INTRODUCTION 

The physical properties of doped and undoped 
graphene are the subject of an intense investigation in the 
view of possible applications of these materials in elec- 
tronic and optical devices. A precise characterization of 
its transport and optical properties is thus a compelling 
issue. In addition to its relevance for technological ap- 
plications, graphene poses several unusual and interest- 
ing theoretical problems. The low energy properties are 
dominated by the Dirac-like excitations at the so-called 
K and K' points of the Brillouin zone, where the tight- 
binding electronic dispersion of the graphene honeycomb 
lattice can be approximated as £k = ±/iUF|k|. Here k 
is the relative momentum with respect to the K (K 7 ) 
point. 1 An additional ingredient is also the chiral char- 
acter of the bands, which gives rise to the well-known 
absence of backward scattering in transport^ Note that, 
although obviously related in graphene, these two is- 
sues are formally distinct and are associated with dif- 
ferent phenomenologies. A Dirac-like behavior can be 
induced for instance in nodal d-wave superconductors, as 
cuprates, without any relation to chirality^^ On the 
other hand a chiral structure is still present even in bi- 
layer graphene, where the low energy bands acquire a 
parabolic character^ 

Transport properties represent a particular delicate 
issue in graphene. As in any other system, also for 
graphene the full quantum treatment of transport pro- 
cesses would require the explicit evaluation of the cur- 
rent vertex corrections for the conduct ivityi^ Indeed, 
only going beyond the so-called bare-bubble approxima- 
tion one can recover the distinction between the trans- 
port relaxation time r tr and the "quasi-particle" one r qp 
(relevant for instance for photoemission measurements). 
However, the chiral nature of excitations in graphene in- 
creases considerably the complexity of this approach. For 
this reason, transport properties are often discussed in 
doped graphene (i.e. when the chemical potential is far 



from the Dirac point) within the framework of the semi- 
classical Boltzmann theory^J&iiil&i&Mii&I&II In this 
context, rtr differs from r qp for a weighted average of the 
scattering probability with the angular factor (1 — cos 8) 
(where is the angle between incoming and outcoming 
scattering electrons). As usual, this leads to the suppres- 
sion of forward scattering processes! 18 ' 19 

Even though a general expectation holds that Boltz- 
mann theory should be valid for doped graphene, the ap- 
plication of this approach to a chiral system as graphene 
is far from being trivial. Indeed, due to the multiband 
chiral structure, the velocity operator Vk = h~ 1 dHi i /dk 
in graphene does not commute with the Hamiltonian 
Hk = hvp\i-er itself, where k = (k x , k y ) and er = (a x , a y ) 
is the vector of the Pauli matrices. As it has been ob- 
served by many authors ) 20 ! 21 ' 22 this fact poses several 
doubts on the applicability of conventional Boltzmann 
theory. In Boltzmann theory indeed one assumes that 
the equilibrium distribution function /(ek) in the pres- 
ence of the external electric field E can be described 
in terms of the one in the absence of external fields as 
/ E (k) w /°(e k - eT tr v k • E), where v k = /rMek/dk 
and e is the electron charge. 23 An important underly- 
ing assumption here is that the energy eigenvalue Ck is 
a good quantum number as well as the shifted quantity 
6k~ er tr Vk-E. This would imply that the energy Hamilto- 
nian operator and the velocity operator Vk commute 
so that they can be diagonalized simultaneously. As dis- 
cussed above, however, this condition is not fulfilled in 
chiral graphene. To overcome this problem alternative 
approaches based on quantum and/or quasiclassical ki- 
netic equations have been employed, where distribution 
functions and the corresponding density operators are 
defined in a chiral matrix spacei 20 ' 21 ' 22 

A second potential limit in the applicability of Boltz- 
mann theory concerns the origin of the angular factor 
(1 — cos6>) in the expression of the transport scatter- 
ing time. Indeed, in conventional systems it originates 
from the momentum dependence of the current operator 
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= (e/h)dH^/dk, which points along the k direction 
in the isotropic case, j(k) oc k/|k|. It is precisely such 
directional dependence which gives rise to the angular 
factor 1- v k - v k < /|v|| « 1 - k • k'/|k 2 | = 1 - cos6> in the 
transport properties^ Things are drastically different 
in the case of graphene where we have 



J'(k) 



e dHk 

h~dV 



evpcr. 



The important feature to be stressed is that the bare 
vertex current operator is here k- independent, i.e. it 
does not depend on the direction of the momentum k = 
(k x , k v ). Thus the possible relevance of the angular factor 
(1 — cos 9) suppressing forward scattering in the transport 
properties should be better assessed. This is particularly 
important when the scattering processes involve a mo- 
mentum dependence more complex that the short-range 
impurity scattering discussed in Refs. l24U25ll26l . A typi- 
cal example is provided by the scattering of electrons by 
acoustic phonons, which plays a major role in controlling 
the temperature dependence of the DC conductivity. 27 ! 28 

In order to assess these open issues we provide in this 
paper an explicit derivation of the DC conductivity in 
doped graphene using a fully quantum approach for the 
electron-phonon scattering based on the Kubo formula. 
Indeed, in non-chiral systems this procedure is known to 
give a quantum derivation of the Boltzmann theory^ 8 - In 
the case of graphene this program can be fulfilled only re- 
taining the full chiral structure in the explicit calculation 
of the current vertex corrections. In the limit when the 
chemical potential is larger than the quasiparticle scatter- 
ing time an analytical solution can be derived. We show 
that, in spite of the above argumentations, at least in 
this regime the results of the Boltzmann theory are still 
valid, in the sense that the forward-scattering suppres- 
sion in the transport relaxation rate is still operative, as 
resulting from the additional angular factor (1 — cos 9) in 
T tr . This result follows from the fact that the dependence 
on the momentum angle of the vertex function is strictly 
replaced by a corresponding angular dependence in the 
pseudospin space, encoded in the Pauli matrix structure. 
We investigate the applicability of the Boltzmann results 
also in the presence of a weak inequivalence between the 
two carbon sublattices, which, close to the Dirac point, 
leads to a mixing of the chiral eigenstates. 29 ' 30 We show 
that even in this case the Boltzmann results are recov- 
ered. The present quantum approach provides in addi- 
tion a basis for the future calculation of the quantum 
vertex corrections at finite frequency and in the Dirac 
limit, which can be both relevant for a direct comparison 
with existing experimental data i 27 i 28 ' 31 

The structure of the paper is the following. In Sec. 
II we introduce the Hilbert space we are working with 
and the basic formalism for the electron-phonon interac- 
tion, while in Sec. Ill we evaluate the self-energy due 
to the electron-phonon interaction in graphene. In Sec. 
IV we implement the calculation of the vertex correc- 
tions for the DC conductivity, and we derive the explicit 



expression for the transport scattering rate. In Sec. V 
we report the results in the presence of an additional 
weak inequivalence between the two sublattices, which 
breaks the chirality preserving the Boltzmann approach. 
In Sec. VI we discuss some possible observable effects of 
our analysis in relation to the doping dependence of the 
high-temperature linear slope of the resistivity, and in 
Sec. VII we summarize our conclusions. In the Appen- 
dices A and B we include some details on the calculation 
of the Eliashberg functions and on the gapped case, re- 
spectively. 



II. THE MODEL 

In order to point out the relevance of the angular trans- 
port factor (1 — cos#) in chiral doped graphene, we con- 
sider in this paper the effects of electron scattering with 
acoustic phonons. 

Let us start by introducing the general electron-phonon 
Hamiltonian in terms of the usual orbital spinor ip^ — 

( c ik A' c ic b)' wnere A , B represent the creation oper- 
ator of one electron with momentum k on the sublattices 
A and B respectively. In this basis we can write: 

k q,l/ 
+ N s X! ^k+qSk.k+q.^k (flq,!/ + O- 



q,/' 



k,q,f 



Here is the non-interacting tight-binding electron 
Hamiltonian, 



Hi 



( o /*(k) 

V /(k) 



(1) 



where /*(k) is the Fourier-transform of the tight-binding 
model on the honeycomb lattice, a q v is the creation op- 
erator of one phonon with momentum q in the v branch, 
cij qjl/ is the corresponding frequency and gk.k+q.^ is the 
electron-phonon matrix element which presents in gen- 
eral a non trivial matrix structure. N s — 2 takes into 
account here the spin degeneracy which will play no role 
in the following. 

In this paper we shall work in the original basis of elec- 
tron operators for the A and B sublattices, as described 
in Eq. ([1]). This choice results to be more convenient 
in the evaluation of the transport properties from a lin- 
ear response theory based on the Kubo's formula. As an 
alternative approach the eigenvector basis, in which the 
Hamiltonian ([1]) is diagonal, was employed in Refs i 20 ' 21 i 22 
since it makes easier the implementation of a quantum 
or quasiclassical extension of the kinetic equations. 

A compelling treatment of the electron-phonon inter- 
action including the full phonon spectrum for all q's is a 
formidable task due to the complex matrix structure of 
the elements gk.k+q.^- Fortunately, since in weakly doped 
graphene the electron momenta are close to the Dirac 



FIG. 1: Diagrammatic representation of the electron-phonon 
interaction. Straight lines represent incoming and outcoming 
electrons, the wavy line is the phonon propagator and the 
filled circles are the el-ph matrix elements. 



FIG. 2: Sketch of the cylindrical coordinates used in the pa- 
per. In this coordinate basis we can write k = fe(cos (j>, sin (f>). 



point, the relevant exchanged phonons are mainly located 
either close to the K, K' edge zone or at the T (q = 0) 
point. Detailed analysis of the electron-phonon effects 
for general acoustic and optical modes can be found in 
Refs. [UlsilaSlliiiiSlsiii]. We concentrate here on 
the electron-phonon scattering with the acoustic modes 
with q w 0, which are dominant in the DC transport 
properties at low temperatures. Several simplifications 
can be employed in this case. First of all, the phonon 



dispersion can be simply linearized, u q 



hv s \q\. In 



addition for q — > the charge modulation wavelength 
A = l/|q| is much larger than the interatomic distance a. 
In this regime the two carbon atoms of the unitary cell are 
essentially indistinguishable so that the electron-phonon 
matrix element behaves as <7k,k+ q ~ <?k,k+ q ^ in the chiral 
space. Within the same assumption, we can also neglect 
intervalley scattering and discuss interactions within a 
single Dirac cone^ Finally, for small doping, we can lin- 
earize the electron dispersion close to the Dirac points 
H k = hv-pk ■ cr. We can write thus our effective Hamilto- 
nian as: 

H = N s N k ^2 ^l^Fk • crip k + ^2 7w s |q|a q a q 

k q 

+N S N K Sq^iUq^k (a q + aL q ) , (2) 

k,q 

where we made use of the relation <?k,k+ q — <? q for q — > 0. 

The resulting electron-phonon coupling is usually ex- 
pressed in terms of the kernel Wk-k'(w — w'), which is 
diagrammatically depicted in Fig. [T]and which represents 
the effective retarded interaction between two electrons 
with momenta k, k' and energies u>, u>', which exchange 
momentum k — k' and energy u> — u>'. From @ we have 



Iflk-k'l -°k-k'(w 



(3) 



where -D q (Sl) = 2aj q /[tj q — 57] is the phonon propagator. 

Note that, since g q oc / in the limit q — > 0, Eq. ([3]) does 
not present a matrix structure, simplifying notably the 
calculations. It is also convenient to express the effective 
retarded interaction is term of the Eliashberg function 



W k -v(w-u') = J dn 
where a 2 F(k - k', n) = \g k -k'\ 



2na 2 F(k-k',n) 



n 2 - (a 
2 8(n - w k _ k 



A2 



UJ>) 



(4) 



In this paper we shall focus on the case of intrinsi- 
cally doped graphene, where the chemical potential 
is much larger than both the allowed exchanged phonon 
energies u; max and the quasiparticle scattering rate r qp . 
Note that in the case of acoustic phonons, the highest 
exchanged phonon energy o> max is given by the Bloch- 
Grimeisen energy scale ma x = 2hv s kp, so that the con- 
straint \fj,\ 3> w max implies vp 3> v Sl which is always 
fulfilled in graphene. The condition T qp <C \fi\, on the 
other hand, is doping-dependent and it is usually ful- 
filled in doped graphene. In this situation, since elec- 
trons are scattered only within a narrow energy window 
±w max around the Fermi level, we can put the electron 
momenta appearing in the Eliashberg function on the 
Fermi surface, so that it depends only on the relative an- 
gle. Writing k = fc(cos </>, sin 0) (see Fig. [2]), we have thus 
k ps k' « kp, and we can write |q| = 2kp sm[((j> — <t>')/2] 
and a 2 F(k - k', Q) = a 2 F(0-0', Q). Moreover, the con- 
dition | /i | ^> o; max , r qp allows us also to restrict ourselves 
to a single electron (hole) cone for fi > (/i < 0), with a 
significant simplification of the calculations [see Eq. (fTS|) 
below] . 

Before discussing the one-particle self-energy, let us 
briefly summarize the properties of the non interacting 
system whose Green's function per spin and valley, in the 
Matsubara space, reads: 



G°(k,iuj n ) 



(ihoJn + n)I — Hvpk ■ cr 
(ihuj n + n)I + hvpk ■ cr 
[{ifvjJn + fi)] 2 — (hvFk) 2 



(5) 



We can expand the Green's function in the Pauli matrix 
basis, 



6°(k,tw„ 



i— I,x,y 



(6) 



It is easy to see, from Eqs. ([5])-([6]), that the diagonal part 
Gj(k, iu> n ) depends only on the e = hv^k (k — |k|) while 
the off-diagonal components depend also on the angle <f>. 
In particular, we can write explicitly 

G°j(k,tuj n ) = G°( Wl ), 
G°(k,ioj n ) = G _(e,ioj n ) cos <j>, 
G° y (k,iuj n ) = G°_(e,iuj n ) sin (j), 
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FIG. 3: Diagrammatic representation of the electron-phonon 
self-energy. Graphical elements as in Fig. [1] 



where 



iflUJr, 



ifkUr, 



III. ONE-PARTICLE SELF-ENERGY 

Let us now consider the one-particle self-energy for the 
electron-phonon interaction. As mentioned, we are in- 
terested here on the limit 3> to max . In this regime 
Migdal's theorem^ assures the validity of the standard 
mean-field like Eliashberg's theory4i The corresponding 
self-energy is thus diagrammatically depicted in Fig. [3J 
and it can be written in Matsubara space as: 

£(k,n) = T Wk-k'{n — m)G(k', m). (7) 



where G(k, n) — G(k, iu> n ), S(k, n) = 

W k -k'(«- - m) = Wk-k'(«W„ - iu> m ). 

From Dyson's equation we can write: 



£(k, iu> n ) and 



G(k,n) 



1 



(ihio n + — /wpk • a — £(k, n) 



Since Eq. ([7]) is a convolution both in momentum and 
frequency space, it is easy to show that the matrix self- 
energy admits the analogous decomposition © of the 
Green's function ([B]), namely 

E(k,n) = E «( k > n )^' 

i— I,x,y 

where E/(k, n) = Sj(e,n), S a (k, n) = £ a (e, 0, n) = 
S off (e,n)cos0, £j,(k, n) = £ a (e, 0, n) = £ ff (e, n) sin </>, 
where the label "off characterizes the o/f-diagonal ele- 
ments of the self-energy. 

Using the Dyson's equation, we can write once 
more Gj(k, n) = G+(fc,n), G x (k,n) = G_ (e, n) cos <f>, 
G y (k,n) = G_(e, n) sin0, with 



G±(e,n) 



1 
2 

±- 



1 



j?Kj„ + /i - e - S/(e, n) - S D ff (e, n) 
1 

iHuj n + + e - S/(e, n) + £ off (e, rt) 



(.8) 



The calculation of Eq. ([7]) is straightforward in the 
case of intrinsically doped graphene. As discussed above, 



the relevant electron momenta are here restricted on the 
Fermi surface, k « kp (e = /i), so that the self-energy 
(J7J) depends only on the angular part of the k vector, 
Ej(e, <f>, n) « Y,i((/) 1 n). We can also split the sum over 
k' in Eq. ([Jj) in its energy and angular degrees of free- 
dom, J2k — y J-n dfi/ZTrf^ deN(e), where W is the 
electron bandwidth, and N(e) — eV/2irh 2 Vp is the den- 
sity of states per spin and per valley, and V is the unit-cell 
volume. Thus we have 



^j(<j),n) = Ty~] / — W<f,-<f,i(n - m)xi((/)')Gi,i oc 



(m), 
(9) 



where \i = h Xx{4>) = cos 4>, Xy{4>) = sm<f>, 
and Gi ! i oc (m) is the local (k-averaged) Green's func- 
tion: Gi : i oc (m) = G +! i oc (m) = / N(e)deG+(e,m), 
G Xy y t \ oc {m) = G_,ioc(m) = / N(e)deG-(e,m). 

It is useful at this stage to introduce the basis of the 
two-dimensional spherical harmonics 

^ a (4>) = e* a t a = 0,±l,±2,..., 

so that we can decompose any generic angle-dependent 
function S(<fi) on this basis, 



where 



S a 



d6 
2^ 



(10) 



(ii) 



Using the definition (01 of the electron-phonon kernel, 
we can write the diagonal and off-diagonal components 
of the self-energy in Eq. ([9]) as 



m J 

E off (n) = T^J 



dCl 



dil 



2n a 2 F (n) 
q 2 + (w n - Lu rn y 

2Qa 2 Fi(fl) 

q 2 + (w„ - Lu m y 



r G+, loc (mX,12) 



r G_, loc (mXl3) 



where a 2 F a are the projections of the Eliashberg function 
a 2 F(4>) on the spherical harmonics ip a (4>), according to 
the decomposition (fT0|) above. 

Eqs. (P~2]) - (P~5|) can be easily analytically continued 
on the real-frequency axis using standard techniques. 
Through the imaginary part of the self-energy, we can 
thus define a diagonal and an off-diagonal scattering rate, 
r /(off) = -lim w ^ ImE /(off) (w + zO + ), where 

Ti = -2nN(n) J dna 2 Fo{n)[n(f3Q) + f(/3n)} 

xIm[G + j oc (f! + i + )] , (14) 

and where 

r off = -2irN(fi) J dfia 2 fi(fi)[n(/?Q) 

xlm [G- Aoc {Q + iO + )] . (15) 
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Here n(x) — l/[e x — 1] and f(x) — l/[e x + l] are the Bose 
and Fermi factors, respectively. 

For a practical evaluation of the self-energy terms, we 
can make use once more of the fact that, due to the low- 
energy phonon-mediated scattering, the relevant electron 
energy are restricted to the Fermi level. In this case 
we can write deN(e) ~ N([i) de and, assuming 
fx > 0, it is easy to see that only the upper Dirac cone 
[first term in Eq. (JSj) ] is relevant. We have in particular 
G+, loc {n+i0+) « G_j oc (Sl+iO+) = -wriV(»/2, so that, 
just as in common metals, the real-part of the self-energy 
vanishes, and we have Ti = 2Kq, T q = 2K\, where 



irN(n) 



d0a 2 F a (0)[n(l3&) + f(P&)}- (16) 



For the states at the Fermi energy, which involve only 
the upper band (/i > 0), we can define thus a total quasi- 
particle scattering rate as 

rqp-rj+Toff = 2tf„ + 2Ki , (17) 

which, using Eq. ©, gives the dressed Green's functions: 



G+(c,uj) w G-(e,ui) w ~ 



2 huj + fi — e + iT qp 
From Eq. |l4])-((T7j) we obtain finally: 



(18) 



i=0A J 

= nN(ri J ^-g 2 e (l + cos 9)[n((30 e ) + f(f3 0e )], 

(19) 

where g$ = I\q\ = I2kp sin(#/2) and 0g = ?iu s |q| = 
2hv s kp sin(6>/2). Using the explicit expression for the 
K a coefficients (see Appendix IA"|) . we obtain standard 
results with a r qp = (7r 2 iV(^)//Sw max z; s )T 2 in the Bloch- 
Griineisen regime (T < ma x) and T qp = (irN(fi)I /hv s )T 
in the high-temperature T > ma x regime. 

It is important to note that the total quasi-particle 
scattering rate arises from both the diagonal and off- 
diagonal contributions of the self-energy. Indeed, the ex- 
plicit angular dependence shown in the second line of 
Eq. pop points out that the off-diagonal terms, giv- 
ing rise to the cos#, are fundamental in order to re- 
cover the usual (1 + cos 9) factor, which accounts for the 
well-known absence of backscattering in graphene due to 
chirality^i^^ 



IV. DC CONDUCTIVITY 

We implement now a full quantum treatment to eval- 
uate the DC conductivity a in the presence of vertex 




FIG. 4: Diagrammatic representation of the the current- 
current response function II. The small empty triangle on 
the left represents the bare current vertex j(k) while the big 
triangle on the right is the renormalized current vertex func- 
tion J(k; n,n + m). 



current renormalization. We consider thus the current- 
current response function^ 



n(m) 



N S N K 
V 



T]TTr[j(k)G(k,n) 



x J(k; n, n + m)G(k, n + m) 



(20) 



where N s , Nk — 2 are respectively the spin and valley 
degeneracies, j(k) = ev-pa x is the bare current operator 
along the x-axis, and J(k; n, n + rn) is the fully renormal- 
ized vertex current. The diagrammatic representation of 
the current-current response function is shown in Fig. 
|U The DC conductivity will be obtained as the limit 
a = — lim w ^oI m II(Lj + i0 + )/uj after the analytical con- 
tinuation of the current-current response function II on 
the real frequency axis. 

It is useful to remark here the importance in transport 
properties of the current renormalization processes which 
account for the backflow cloud associated to the current 
of the quasiparticles. As discussed in the introduction, 
one of the main effects of such current renormalization 
in normal metals is to give rise to the angular factor 
(1 — cos 9) which differentiates the transport scattering 
rate r tr from the quasi-particle one r qp . We can under- 
stand this result by noting that, in normal (isotropic) 
metals, the bare and the renormalized currents, j and 
J respectively, are both proportional to the velocity, i.e. 
they point in the direction of the momentum k. Thus, 
one can write J — Ak, where the vertex function L can 
be computed for instance in the ladder approximation. 
The self-consistent solution gives thus£ 



1 



kk' 



(21) 



where k, k' are the momenta of the scattered electrons, 
with modulus equal to fcp but different directions, and 
(■■■)e indicates the angular integration. This leads to 
the usual additional factor (1 — k-k'/fc 2 ) = (1 — cos0) 
in the angular average of the transport scattering time, 
which reproduces by means of a full quantum treatment 
the well-known semiclassical Boltzmann result. 

While deriving the result (|2"Tj) . a crucial ingredient is 
the proportionality between J and the momentum k. In 
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graphene, where the energy-momentum dispersion is lin- 
ear, such a relation clearly does not hold. However, as 
it is evident already at the level of the bare current j, 
the matrix structure plays the analogous role of the mo- 
mentum dependence in ordinary metals, so that j x for 
example is a matrix proportional to a x . As far as the 
renormalized current is concerned, the analogous of the 
momentum dependence in the ordinary metals becomes 
now a decomposition of J in the Pauli-matrices compo- 
nents, by means of dimensionless function A defined by 
the relation J(k; n,n + m) = eupA(k; n, n + m). As we 
shall see in what follows, in graphene such a matrix struc- 
ture compensates - in a non-trivial way - the lack of mo- 
mentum dependence of the Fermi velocity, and leads once 
more to the Boltzmann result. 

Aiming on focusing on the matricial structure of the 
current function, in the following we shall make use for 
this quantities of the same approximations employed for 
the self-energy. In particular, taking into account that 
the electron-phonon interaction gives rise only to low- 
energy scattering, we can approximate k ~ fcp in the 
vertex function and retain only the angular dependence. 
In addition, one can decompose the vertex function in 
the basis of the Pauli-matrices, 



A(</>; n,n + m) 



Ai(4>;n,n + m)ai, (22) 



and we can expand the Aj(</>) functions in terms of 
the spherical harmonics components as in Eq. (JTUJ): 
Ai(0; n,n + m) — J2 a ( n > n + TO )V'a(^)- Inserting Eq. 
([22)1 in Eq. ([20]). and assuming once more, for doped 
graphene, G+ « G_ (neglect of inter-band scattering), 
we obtain thus the general structure: 

II(m) = ^^T^^in, n + m)b(n, n + m), 



V 



(23) 

where b(n, n + m) = N(fi) J deG(e 1 n)G(e, n + m) (since 
G+ = G_ we have further dropped the index "±") and 
where L tot (n,n + m) = J2i, a cf L?(n, n + m). The cf 
are numerical coefficients which arise from the angular 
average over </>. Since in the bubble (|20|) one has the 
product of two Green's functions, it is easy to realize 
that the angular average will involve at most harmonics 
up to the second order. Indeed, an explicit calculation 
shows that the only non vanishing terms are 



1, 



-ic„ 



(24) 
(25) 



Once the analytical continuation iu) m — > u> + i0 + of the 
bubble (|2"3"|) is performed and the limit id — > is evaluated 
to compute the DC conductivity, one can see^ that only 
the advanced-retarded parts of the vertex function and 
of the 6-function contribute, so that: 



2, ,2 



NfcN s he v 



-6 RA (0)L^(0), 



(26) 



A 



> 



+ 



A 



FIG. 5: Diagrammatic representation of the current vertex 
function. Graphical elements as in Fig. [T] The ladder ap- 
proximation is enforced by the Migdal's theorem valid for 
\n\ » w max ■ 



where 6 RA (e) = lim w _ b(ui+i0 + , uj-i0 + ), andL^X(O) = 
lim w ^ h tot (w + i0 + ,w - i0 + ). From Eq. (fill) one easily 
finds 



6ra(0) 



4r nn : 



which gives the general expression for the DC conductiv- 
ity in the presence of vertex corrections: 



N K N s he 2 vlN{n) 
4VT tr 

where we define the transport scattering rate 

r — 1 qp 
tr "L-(0)- 



(27) 



(28) 



In the absence of vertex corrections (Ijra(0) = 1) one 
recognizes in Eq. ([27]) , as mentioned above in the case of 
ordinary metals, the standard result of the conductivity 
in the bare-bubble approximation ; 24 ! 42 where the scatter- 
ing time for transport coincides with the quasi-particle 
scattering time a = e 2 /i/(27rftr qp ). 

To compute Eq. (|2"6")l we need then to calculate the 
vertex corrections. In the doped graphene regime we 
are interested in, this aim is made easy once more by 
the Migdal's theorem which enforces the validity of the 
mean-field theory. The corresponding vertex function 
A(k; n,n + m) can be evaluated thus within the self- 
consistent ladder approximation (see Fig. [S]), namely^ 

A(k; n, n + m) = a x + T V" W k _ k ' (n — l) 

k',Z 

xG(k', Z)A(k'; /, / + rn)G(k', I + m). (29) 

In the same limit ^> w ma x, we can also employ the 
approximations implemented for the self-energy, namely 
|k| « |k'| « fc F , — N^J^d^^J^de', and 

we can also set A(k F ; I, l+m) — A(<p'; I, l+m) in the right- 
hand side of Eq. ([2^1 . Expanding again the A(</>'; I, l+m) 
in terms of the Pauli basis and of the spherical harmonics, 
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we end up with the following set of equations 
A^(n,n + m) = S afi + c^W a (n - I) 



xb(l,l + m)A^ +a {lJ + m), (30) 



Aj(ra, n + m) 



A"(n, n + m) 



1,3,0 



xb(lJ + m)A^ +a (l,l + m), (31) 
Tj2fjW a (n-l) 

1,3,0 



xb(l,l + m)A? +a (l,l + m), (32) 

where the coefficients cf are defined in Eqs. (|2"4")) - (f2l))) . 
and the only non zero terms of the coefficients df, /" 
are: 



d°j = 2, 

d\ = dZ 1 = -id\ = 1 

fy = -if} = iff 1 = 1, 

f2 



2 



We consider now the analytical continuation of Eqs. 
(f3"0")) - (|32p on the real axis. We are interested on the quan- 
tity A" RA (0). We can apply the usual standard proce- 
dures for the analytical continuation 8 to each element on 
the right side of Eqs. (|3"0"|) - (|3"2"f . and we get, for instance, 



A% A (0) = S afi + 



tt7V(m) 



dn a 2 F a (n) 



x[n(fi(l) + f(/Xi)] 



v 

^ 3 r qp (^) 



3,0 



(33) 



Similar expressions hold true for Aj RA (0), A^ RA (0). 

The quantities A^ RA (fi), r qp (fi) in the right side of 
Eq. (f33|) have a significant fl variation over an elec- 
tronic range of energies, whereas the Eliashberg function 
a 2 F a (Q) limits the energy integration up to the phonon 
scale w m ax- In this energy range we can then approxi- 
mate in the integral of Eq. A^ RA (Q) ~ A^ RA (0), 
r qp (f2) ~ r qp (0) = r qp , and we obtain a simple set of al- 
gebraic relations for the vertex function at zero frequency, 



AWO) 



A^ra(O) 



r, 



'V<--^£cfA£ RA (0), 



qp 



3,0 



qp 



3,0 



qP 3,0 



Despite its apparent complexity, such a system of 
equations admits a simple solution. In particular, by 



exploiting the symmetric/antisymmetric properties for 
a — ► —a of the x,I and y components, respectively, 
one can see that only three independent components are 
not zero, i.e. 



A° RA , A 2 . RA , A} RA , with the relations, 



A 2 



RA 



A~ 2 — i A 2 — — ?'A~ 2 A 1 



= A7,b 



RA- 



The system can be further simplified by noting that it 
can be rewritten in terms of a single self-consistent equa- 



tion for A RA : 



Ar° a = 1 



RAi 



whose solution gives 

AM = 



qp 



r qp - K - 2Kt - K 2 



(34) 



Finally, by expressing the quasi-particle scattering rate 
r qp as a function of the K a terms as in Eq. (j 1 T[) we 
obtain: 



T tr = K*-Ko 



(35) 



Eqs. (|34|). (|35|) are the main result of the present pa- 
per. Their physical insight appears clearly if we express 
the transport scattering rate in terms of the microscopic 
electron-phonon interaction, in analogy with the result 
(fT9|) for the quasiparticle scattering rate: 



Ttr 



K - K 2 

, T/ , f d9 ,1- cos 20 r ,„ . S1 
/ 2^ 9e 2 [ n (fae) + f(p0e)] 



-kN(h) 



dB 
2^ 



g 2 (l - cos 6»)(1 + cost 



x[n(P0 ) + f(P0 9 )}. 



(36) 



Taking into account the expression of the quasi-particlc 
scattering rate r qp in Eq. (fT9|) . we can write Eq. ([36]) as 



(r tr ) e = <r qp (i-cos0)) e . 



(37) 



This analysis shows thus that a fully quantum deriva- 
tion of the DC conductivity yields in doped graphene 
the same result than the standard Boltzmann theory. In 
particular we can see that, although the vertex function 
is k-independent in graphene, the effect of current ver- 
tex corrections is to add, just as in common systems, 
the additional angular factor 1 — cos 8 in the phase space 
probed by the electron-phonon interaction, giving rise 
to a suppression of the forward scattering. The effects 
of such angular factor can be remarkably traced in the 
temperature dependence of the transport scattering rate 
r tr = Kq — K 2 (and hence of the resistivity p = 1/er) com- 
pared with the quasi-particle one r qp = 2(Ko + Ki). An 
analytical derivation of the K a coefficients is reported in 
Appendix [X] and the corresponding temperature depen- 
dence of r tr and r qp is shown in Fig. [5] These results 
are pretty well consistent with the standard Boltzmann 
theor y 18 i 19 In the high temperature limit T ^> max only 
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FIG. 6: (Color online). Temperature behavior of the quasi- 
particle scattering rate F qp and of the transport scattering 
rate I\r for p, — 0.1 eV, corresponding to m ax ~ 4 meV. The 
transport scattering rate F tr can be also expressed as a func- 
tion of resistivity, p = AWxy/Nk N B he 2 VpN(p), as plotted on 
the right-side scale. 



the Kq component survives, with Kq = irN(p)IT /2hv s 
(see Appendix E|, so that r tr ,r qp ~ T. On the other 
hand in the Bloch-Griineisen regime, just as in common 
metals, we recover the usual r tr ~ T A while T qp ~ T 2 . 
This different dependence can be understood by noting 
that the low temperature behavior of T qp in Eq. (| 1 7|) is 
dominated by the leading order ~ T 2 of K$ , K\ , whereas 
the leading orders Kq, K2 in Eq. ([55)) cancel out so that 
the temperature behavior of r tr in this regime stems from 
the higher order ~ T 4 contributions. 



V. BREAKING THE CHIRALITY 

In the previous Section we have evaluate the DC con- 
ductivity in chiral doped graphene by using a quantum 
approach based on the Kubo formula. We have shown 
that, although the vertex function is momentum indepen- 
dent, the Boltzmann results are fully reproduced, even 
for what concerns the presence of the so-called angu- 
lar transport factor 1 — cos 9 which suppresses forward 
scattering. We would like to stress however that the re- 
sulting validity of the Boltzmann theory is by no means 
trivial. Indeed, as we have shown, the robustness of the 
Boltzmann results is due to the chiral structure itself of 
graphene, = hv F k • <x, which translate in the direc- 
tion of the Pauli matrix pseudospin <x the role which in 
common metals is played by the momentum k-direction. 

Such observation rises the question about the valid- 
ity of the Boltzmann result when the chiral structure 
of graphene is affected by external fields. This hap- 
pens for instance in the presence of a weak sublattice 
inequivalence, where an energy potential difference be- 
tween the two sublattices, besides opening a gap at the 
Dirac point ) 34 i 43 gives rise also to a mixing of the chiral 



eigenstates and to a loss of chirality close to the Dirac 
point. Such reduced role of the chirality close to the 
Dirac point has been traced, for instance, in Ref. I29H30I 
in the intensity profile along the constant-energy contours 
probed by angle-resolved photoemissioni 29 ' 30 The valid- 
ity of the Boltzmann results even in this context needs 
thus to be revised. 

To this aim in this section we consider doped graphene 
in the presence of a weak inequivalence between the two 
carbon sublattices, parametrized in terms of a different 
sublattice potential A: 

i?k = hv F k- a + Aa z . (38) 

The energy spectrum is easily obtained from (|38|) . E% = 



yj (hvpk) 2 + A 2 , showing the opening of a gap at the 
Dirac point. 

Let us consider for the moment the case of a non- 
interacting system. The presence of A induces an 
explicit cx a z component in both the Green's func- 
tion and self-energy. Eq. ([6]) is thus generalized 

as G(k,n) = Y,i=i x y z G i(^ n )^i' wnere Gi(k,n) = 
G+{E,n), G x (k,n) = \' oS (E)G-{E,n) cos<f>, G y (k,n) = 
j oS (E)G-(E,n) sin0, G z (k, n) = 7z (E)G-(E, n), with 



G±(E,n) 



1 



± 



1 



ihuj n + fj, — E ifku n + ji + E 

Joft(E) — e/E, j e (E) = A/E. Here we denote e = hv-pk 
and E = y/e 2 + A 2 . 

Since we are interested in the doped graphene, where 
\fi\ ^> w max , and we assume pi > 0, we can as usual 
consider only the upper band, and we can restrict the 
electr onic states on the Fermi surface, for which E as p, 
e ~ yl - A 2 //i 2 , and we have 



Gi(E,n) = ^g(E,n) 



(39) 



where 7/ = 1, j x = j y = 7 off = y/l - A 2 /p 2 , ^ z = A/p 
and where g(E, n) = l/(ifko n + p — E). 
Plugging Eq. (f3"9"| in 0, we can write 



Ej(n) 
E fi(n) 
E z (n) 



^2 Wi(n - m)g loc (m), 
m)g ioc (m), 



where g\ oc (m) — N(p) J dEg(E,m). After the standard 
analytical continuation on the real-axis, we can see as 
usual that the real part of self-energy vanishes and the 
imaginary parts of the different components give rise to 
the corresponding scattering rates 



r, = 2K , 

r ff = 7off2id, 
r 2 = lz 2K Q . 



(40) 
(41) 
(42) 
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Inserting this relations in the matrix expression for the 
Green's function, after some careful derivation described 
in Appendix [Bl we can write a Green's function in the 
presence of interaction in the form of Eq. (fT5)l . with 



r qp = Tj + 7offr o ff + -y z T z 
= (1 + 7^)2^0 + 7^2^ 



A 2 

M 2 



1K X . 



We can write thus 



qp 



A 2 
M 2 



Eq. (|4"3")l shows in a direct way the reduced effect of the 
chirality on the quasi-particle scattering rate in the pres- 
ence of a sublattice inequivalence. For instance we can 
see that the back-scattering 9 — it is now not completely 
suppressed but it is reduced by a factor A 2 /p 2 with re- 
spect to the forward-scattering 6 = 0. For \p\ = A, when 
the chemical potential is on the edge of the gap^ the 
K\ angular dependence disappears completely and the 
quasi-particle scattering is completely isotropic with no 
effect of chirality. 

Let us consider now the DC conductivity. Using the 
same derivation as in Sec. IIV1 we can still write a set 
of equations as (j37))) - (j3"2")) with the same cf, dfj" 1 and 
rescaled quantities K — ► T%gKa> K\ — * (1 + lz)Ki, 
K2 — > 7 2 ff i^2, The presence of cx & z terms in the Green's 
function gives rise in addition to a finite z-component 
A". After lengthly but straightforward calculations we 
can still write the same self-consistent equation for the 
total vertex function as before, with the rescaling of the 
K a coefficients discussed above. The final result is 



ArI(o) 



7 2 ffr qP 



so that 



Itr 



r qP - 7 2 ff ^o - 2(1 + 7l)^i ~ 7 2 ff K2 ' 



T qp - 7 o 2 ff iT - 2(1 + 7 2 )A-! - 7 2 ff K 2 
7 2 ff 

M 2 + 3A 2 4A 2 

■AO 5 ^T-ftl — K 2- 



p 2 - A 2 " p 2 - A 2 



(44) 



Eq. ([44|) looks at a first sight quite dim. However, tak- 
ing into account the explicit angular dependence of the 
Eliashberg functions a 2 F a and of the coefficients K a , we 
can re-arrange Eq. (|44|) in the form: 



(r tr ; 



1-^2 ) 'T„„(l -c(..s«);„. 



(45) 



where the presence of the angular term (1 — cos 9) sug- 
gests again the validity of the Boltzmann result even in 



this case where the chirality of graphene close to the 
Dirac point is affected by the sublattice inequivalence. 
As far as the prefactor is concerned, one could in prin- 
ciple justifies this as well within the Boltzmann frame- 
work by taking into account that the the opening of the 
gap affects also the Fermi velocity of the eigenstates of 
the Hamiltonian (|3"8j) . Thus, following the Boltzmann 
prescription, the velocity in the basis of the eigenstates 
is v k = h- 1 d^J{hv ¥ k) 2 + A 2 /dk = v F e k /yJe{+ A 2 w 
A 2 //i 2 , where e k = hvpk and y/ ej. + A 



/' 



v FX /T 

at the Fermi level. Thus, the additional (1 — A 2 /// 2 ) 
factor, that we obtained in Eq. I|45p associated with ver- 
tex corrections, will appear in the Boltzmann expression 
for the conductivity ctb = e 2 N([i)v 2 /2{Tb) associated to 
the v 2 term instead of the scattering rate, which is simply 
(43) (T B ) fl = (r qp (l-cos0)) 9 . 



VI. DISCUSSION 

In the previous section we have shown that, in the 
limit where the chemical potential is the largest energy 
scale of the system, the Boltzmann theory is still valid 
for graphene even in the presence of a weak sublattice 
inequivalence which gives rise to a gap A. A possible 
outcome of the extension of these calculations beyond 
this "Boltzmann regime" is suggested by the observation 
of the (weak) doping dependence of the resistivity p that 
follows from Eq. ()45[> . As it is shown in Fig. [7^,, both the 
low temperature (T <C m ax) and the high temperature 
(T ^> max ) regimes acquire a doping dependence, here 
parametrized in terms of the ratio r = (A//1) 2 . Partic- 
ularly interesting appears the high-temperature regime 
where, using Eq. (I44|) and (|A4j) . one finds: 



P = 



3A 2 



ttD 2 



p 2 — A 2 ihe 2 VpV 2 p, 



-T, 



(46) 



which implies a significant increase of the linear slope 
as one approaches the Dirac point p — > A, as shown 
in Fig. [7Jd. Observe that, if one would use instead the 
single-bubble result for the DC conductivity, the resistiv- 
ity would be proportional to the quasiparticle scattering 
rate r qp , so that according to Eq. ([4"3")) one would obtain: 



Pba 



1 , A 2 



nD 2 



4he 2 v 2 ,v 2 



F v sPm 



-T, 



(47) 



which has a much weaker doping dependence, as shown 
in Fig. [7)3. Remarkably, an increase of the linear slope as 
p decreases has been observed in recent measurements of 
resistivity in doped graphene samples! 27 ' 28 Indeed, while 
the measured crossover from a power-law to a linear T- 
behavior points toward the electron-phonon scattering as 
the source of the temperature dependence, the fact that 
the slope of the linear term decreases as doping increases, 
and eventually saturates at enough large doping, is un- 
explained within existing theories. Eq. (|4"6")) could sug- 
gest a possible mechanism for this doping dependence, 
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FIG. 7: (color online) : (a) Temperature dependence of the re- 
sistivity for different doping in the presence of a weak sublat- 
tice inequivalence, giving rise to the gap A. Here r = A 2 //i 2 . 
Resistivity has been computed using the transport scattering 
time in Eq. (|44[) . Inset: same quantities on a log-log scale, 
(b) Doping dependence of the high-temperature slope of the 
resistivity in graphene in the presence of a gap ( A = 5 meV) . 
The dashed (red) line is the result (|47[) obtained with the bare 
bubble, while the solid (black) line represents the resistivity 
(|46|) computed with full inclusion of vertex corrections. The 
shaded area at [i > 30 meV represents the range of doping 
investigated in Ref. |28| . while in the grey area around 5 meV 



the validity of Eq. (|46[) fails because 
meV has been computed here for T - 



A < r qp (r qp = 1.4 



300 K). 



even though in the regime of validity of this equation the 
quantitative variations of the slope are expected to be 
much smaller than what experimentally measured. In- 
deed, recent tunneling experiments performed in a simi- 
lar suspended graphene sample have found a gap at the 
Dirac point A ~ 5 meV, 44 . This small value of the gap, 
compared to the relevant values \x ~ 30 — 70 meV of the 
chemical potential in the samples investigated in Ref. [H, 
would lead, in Eq. (j4"rj)) , to a variation of the slope of only 
few percents, as shown in Fig. [TJa. However, one should 
keep in mind that Eq. (f4"rJ)) has been derived in the limit 
r qp -C /i. Thus, only the full solution in the crossover 
region where the Dirac point is approached can discrim- 
inate if the doping dependence reported in Refs. I27ll28l is 
an effect of a small gap opening which goes beyond the 



Boltzmann approach. Analogously, it can be worth also 
to investigate alternative scenarios^ 5 - for the gap open- 
ing that have been suggested by photoemission results in 
epitaxial graphene ; 46 ' 47 where the gap A ~ 100 meV is 
comparable to the chemical potential in the interesting 
doping regime. 



VII. CONCLUSIONS 

In summary, in this paper we outlined the issue of the 
calculation of the current vertex corrections for electron- 
phonon scattering in graphene within the Kubo formal- 
ism. While previous works investigated the role of vertex 
corrections in the presence of impurity scattering ) 24 ' 25 ' 26 
here we address for the first time the issue of vertex 
corrections in graphene when the scattering mechanism 
arises from electron-phonon interactions, which have a 
non-trivial momentum dependence. We analyzed in par- 
ticular the case of doped graphene, when [i is the largest 
energy scale of the system (i.e. |/z| 3> w ma x, r qp ). In this 
regime the calculations can be performed explicitly, lead- 
ing to an analytical derivation of the transport scattering 
rate which appears in the DC conductivity. Remarkably, 
we found that despite the lack of direction-dependence of 
the quasiparticle velocity in graphene, the matrix struc- 
ture associated to the current vertex plays a similar role 
than the momentum dependence of the renormalized ver- 
tex in ordinary metals, leading to a confirmation of the 
Boltzmann approach. Such a result is confirmed also in 
the presence of a sublattice inequivalence, which leads to 
a gapped Dirac spectrum and introduces a doping depen- 
dence of the resistivity that can be a promising candidate 
to explain existing experimental data in doped graphene 
samples. With respect to the formalism based on the 
quantum kinetic equations ; 20 ' 21 ' 22 which is also aimed 
to investigate the applicability of the Boltzmann result, 
the present approach based on the calculation of vertex 
corrections for the conductivity bubble has the main ad- 
vantage of being suitable of a direct extension at finite 
frequency. Such an extension will lead to a full quantum 
treatment of the conductivity in graphene, and will be 
the subject of a future work. 
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APPENDIX A: HARMONIC COMPONENTS OF 
THE ELIASHBERG FUNCTIONS AND K a 
COEFFICIENTS 

In this Appendix we provide an explicit expression 
for the harmonic components of Eliashberg functions 
a 2 F a (Q), and we calculate the temperature dependence 



11 



of the corresponding functions K a (T). We recall the def- 
inition of the momentum-dependent Eliashberg function, 

a 2 F{k - k', Q) = | 5k _ k , \ 2 5(Cl - Uk-v), 

where w q w 7jv B |q|- Following Ref. |l9|, we notice that 
|ffq| 2 oc |q|, so that we can write |g q | 2 = J|q|, where 
/ = hD 2 /2Vp m v s , D is the deformation potential, V is 
the volume of the two-dimensional unit cell, and p m is the 
graphene mass density. Typical values are D ~ 19 eV, 
p m = 7.6xl0~ 8 gr/cm 2 , v s = 2xl0 6 cm/s, v F = 10 6 m/s. 
Restricting the electron momenta on the Fermi surface, 
|k| |k'| as fcp, we can write |q| = 2fcp sin[(0 — <f>')/2], 
and the Eliashberg function will result to depend only on 
the exchanged angle 9 = <f> — </>': 

a 2 F(0, fi) = 2Ik F sm(9/2)6[n - w max sin(0/2)], (Al) 

where we remind w max = 2hv s k F . 

The angular components a 2 F a will be evaluated sim- 
ply projecting the Eliashberg function (|A1|) on the spher- 
ical harmonics ip a (4>) — e m< ^. We have thus: 



a 2 F a (0) = 2Ik F [ —sin{8/2)e %ae 5[n-uj n 
Jo 2tt 



send the upper limit of integration to infinity, f3uj n 
oo. We obtain thus: 



K a = 



N(p)I T 2 



a i 



a 3 - 



1 -Aa 2 T 2 

2 0max 



where a n = J °° dyy™[n(y) + /(y)], and in particular, a\ = 
7r 2 /4 and 03 = ir 4 /8. In the opposite high-temperature 
limit T 3> 0max we can instead compute the integrand 
function in Eq. (|A3j) as y — * 0. In this case, due to the 
rapid oscillations of the cosine term, all the harmonics 
vanish except the a = one, which displays a linear T 
dependence: 



(A4) 



APPENDIX B: TOTAL QUASI-PARTICLE 
SCATTERING RATE IN THE PRESENCE OF 
SUBLATTICE INEQUIVALENCE 



2Ifa 



71 dO 



sm0e i2ae 6\n 



In this Appendix we derive the effective total quasi- 
sin(0/2)] particle scattering rate in graphene in the presence of 
sublattice inequivalence. 
gm 0i Let us start from the scattering rates in Eqs. (jit})) - (|4"2"|) , 

which we can write in the matricial form: 



The 5-function has two solutions for 9 6 [0 : n] 
one for 6 = yn and one for 6 = ir — yn, where yn = 
arcsin(r2/w max ). We can write thus: 



a 2 F a (n) = 21 k F f —sm9e aae - 



1 



x [5(9 - y n ) + , 
AIk F fl9(uj ma 



!^ max cos(6')| 
- 7T + yn)} 



n) 



1 



n 



1 r 



e i2ayn _i_ e i2a(n-ya) 



21 0cos(2ayn)9(uj n 



fi) 



(A2) 



where we made use of sin(yn) = sin(7r — yn) = fl/u! max , 
I cos(yn)| = I cos(tt - y n )| = y/l- (Q/uj max ) 2 . 

Inserting Eq. (|A2|) into the Eq. p^|) which defines the 
if Q function, we can obtain their temperature depen- 
dence. By rescaling the integration variable to y — (30 
we have: 



N(p)I T 2 
x cos I 2aarcsin 



dy 



y 



(3uj„ 



V 1 ~ (y//3^max) 2 

[n(y) + f(y)} . (A3) 



In the low-temperature limit T <C m ax we can retain the 
leading terms in powers of yj /30 max in the integrand and 



T = Yil + r off [cos <pa x + sin0<7 y ] + T z & z . (Bl) 
Considering the Hamiltonian in Eq. (j38|) , we have then 

G- l {e,u) = {hiu + iT I )i+{A-ir z )a z 

— (e — «r o ff ) [cos 4>cr x + sin <\>b y ) . 

The spectral function, whose width determines the total 
effective quasi-particle scattering rate, is associated with 
the Gj term which results: 



Gx(e,a,) = 



hw + p, + iTj 



(fiw + m + ir/) 2 - (e - ir off ) 2 - (A - i 
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In order to identify the total effective quasi-particle 
scattering rate, we expand now the denominator and 
write, in particular,: 

(e - zr off ) 2 + (A - iT z ) 2 



a 2 - r 2 ff 



r 2 -2i( e r off + Ar z ). 



Since e 2 + A 2 w p and T g,T z <^ \p\, we can neglect 
T 2 ff — T 2 with respect to E 2 = e 2 + A 2 and, at the same 
order, we can write: 



E 2 - 2i{eT oS + AT Z ) 



_ 2 



A 2 „ A 

— rr ff h — r z 
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Plugging this result in (|B2[) , we can now split the re- so that at the Fermi level in the upper band we obtain 
suiting Green's function in the usual two contributions 
from the upper and lower band, 

Gi(E,u) = lj2 

with 



2^ + ^ + iT lT E±iT^ r qp = r / + ,/i-^r off + ^r 2 



r ± = r /± (Ji-^ ratt + ^r, V 
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